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In conventional approach to the chameleon mechanism, by assuming a static and spherically 
symmetric solutions in which matter density and chameleon field are given by p = p(r) and (j> = (t>{i~), 
it has been shown that mass of chameleon field is matter density-dependent. In regions of high matter 
density such as earth, chameleon field is massive, in solar system it is low and in cosmological scales 
it is very low. In this article we revisit the mechanism in cosmological scales by assuming a redshift 
dependence of the matter density and chameleon field, i.e. p = p(z), (f> — (j>(z). To support our 
analysis, we best fit the model parameters with the observational data. The result shows that in 
cosmological scales, the mass of chameleon field increases with the redshift, i.e. more massive in 
higher redshifts. We also find that in both cases of power-law and exponential potential function, 
the current universe acceleration can be explained by the low mass chameleon field. In comparison 
■ with the high redshift observational data, we also find that the model with power-law potential 

' function is in better agreement with the observational data. 
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I. INTRODUCTION 



o 
a" 

^ , Recently, the observations of high redshift type la supernovae and the surveys of clusters of galaxies Q-Q reveal 
(3JT)' cosmic acceleration. Further, observations of Cosmic Microwave Background (CMB) anisotropies indicates that the 
universe is fiat and the total energy density is very close to the critical one p| . 

These data properly complete each other and designate that the so called dark energy (DE) is the dominant 
| component of the present universe, occupies about %73 of the cosmic energy , while dark matter (DM) occupies %23, 
. and the usual baryonic matter takes about %4. There are prominent candidates for DE such as the cosmological 
Q^ ' constant [1, 0] , a dynamically evolving scalar field ( like quintessence) @, Q or phantom field ( with negative energy) 
C , [Io| that explain cosmic accelerating expansion. Alternatively, the universe acceleration can also be elucidated through 
*n ■ modified gravity brane cosmology and so on [H|]-[l8j]. The DE can track the evolution of the background matter 
\Q \ in the early stage, and only recently, it becomes dominant . Thus, its current condition is nearly independent of the 
O ■ initial conditions [19j|-[23|]. 

' On the other hand, to explain the early and late time acceleration of the universe, it is most often the case that 
such fields interact with matter; directly through matter Lagrangian coupling, or indirectly by coupling to the Ricci 
scalar as emerged from quantum loop corrections [24|-[28|. If the scalar field self-interactions are negligible, then the 
experimental bounds on such a field are very strong; requiring it to either couple to matter much more weakly than 
gravity does, or to be very heavy [2^]-[32j. Unfortunately, such fields are usually very light and its coupling to matter 
should be tuned to extremely small values in order not to be conflict with the Equivalence Principal [33[ . 

An attempt to overcome the problem with light scalar fields has been suggested in chameleon cosmology [34|- 
[3fj| . In the proposed model, a scalar field couples to matter with gravitational strength, in harmony with general 
expectations from string theory whilst at the same time remaining very light on cosmological scales. The light field 
on cosmological scales is permitted to couple to matter much more strongly than gravity does, and yet still satisfies 
the current experimental constraints. The cosmological value of such a field evolves over Hubble time-scales and could 
potentially cause the late-time acceleration of the Universe |37j • In this approach the mass of the scalar field depends 
on the local background matter density. While the idea of a density-dependent mass term is not new [38|-[43j|, in the 
work presented in [3(| [I?) the scalar field couples directly to matter with gravitational strength. In this paper, we 
stat with the assumption that the matter energy density is redshift dependent, p = p(z). We therefore find that the 
chameleon scalar field has a greater contribution to the dynamic of the universe in higher redshifts. 

The manuscript is organized as follows. In section two we describe the chameleon mechanism. We assume two 
forms of power law and exponential potential functions and derive the required equations. Section three is devoted 
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to constrain the model parameters with the observational data using x 2 statistical method. In section four with 
numerical calculations and best fitted model parameters, we revisit the chameleon mechanism and test the model 
against recent findings. Summary and remarks are given in section five. 



II. THE MODEL 



We begin with the action of chameleon gravity given by, 

/M 2 1 
l-^K - ~W + V(<f>)]^dx 4 

+ J £ m (ip (l \g$)dx\ (2.1) 

where the matter fields ip^ are coupled to scalar field <p by the definition gfy = e 2 ^ i ^^ Mpl g^. The (3i are dimensionless 
coupling constants, one for each matter species. In the following, we assume a single matter energy density component 
p m with coupling j3 [36j . The variation of action (|2.1j) with respect to the metric tensor components in a spatially flat 
FRW cosmology yields the field equations, 

3H 2 M 2 pl = Pm e^ + ±^ + V (</>), (2.2) 



(2H + m 2 )M 2 pl = - 7 p m e*V - I</, 2 + V{4>). (2.3) 

In deriving the field equations we assumed a perfect fluid for matter field with p m — jp m - Variation of the action 
(|2.ip with respect to scalar field <j> provides the wave equation for chameleon field as 

4> + ZH4> = -V' -^- Pm e^, (2.4) 
M p i 

where prime indicated differentiation with respect to 4>. From equations (|2.2I) . (|2.3|) and (|2.4[) . one can easily arrive 
at the relation 

Prn + 3Hp m (l+"f) = -3j-^—p m 4>- (2.5) 

M p l 

Integrating the above equation yields 

Pm = > (2-6) 

e 15 ^" a 3 ( 1+ ^) 

where A as a constant of integration. In the following we discuss two power-law and exponential forms for the potential 
in the model. 
Case 1: 

The runaway inverse power-law potential, also called the Ratra-Peebles potential, in chameleon cosmology is often 
given as 

M i+a 

where M is a constant with the dimension of mass and a is a positive constant. These kind of potentials are usually 
seen in quintessence models. The chameleon effective potential is then defined by, 

V eff ^) = V(cj ) )+p m e^, (2.8) 

the sum of potential (|2.7[) and its coupling to the matter density. Using equation (|2.6[) , the effective potential can be 
rewritten as 

]\/f 4 + a (1-37)3 j. 

Veff {ct>) = ^L^+Ae^^^a- 3 ^. (2.9) 



3 



One can also easily find the mass associated with the field </> from: 

™ 2 = W VeU{4)) ' (2 ' 10) 

For matter dominated universe, 7 = 0, and positive j3 in the second term of V e // (</>), the effective potential monoton- 
ically decreases to a minimum at a finite field value <f> = (f> m ini where ■jxVef f\(t>=^ mi „ = 0, and m = m m i n . We then 
lead to the following relation 

= p (2.11) 
where A = (1 - 37)-^ and d = "^ 4+ y 1+T) . From equation ([2TTT]) we find 

<j>min = \w(\(j)?^), (2.12) 
where W is called Lambert or product-log function. We can also find minimum mass from 

m 2 min = a(a + + 4AV*"*» a" 3 ^. (2.13) 

rmin 

Using equation (|2.11|) we can rewrite equation f)2.13p as 

mL„ = aM 4 +«(^±i + -A T ). 

mm r'min 



f4+c 3(1 + 7) 



From a = ^7 we find <i = ( 1+z y?(i+-y) where e?o = — — — is a positive constant. Thus, </> m , ra can be rewritten as 

-bmin = ^(A[^(l + z)-^]Th) (2.15) 



A v A 



Case2: 

The exponential potential given by 



V(cj>) = AfV^, (2.16) 



is also seen in many cosmological models. Similar to the previous case, we find the effective potential, V e ff (</>), the 
minimum mass m m in, and the minimum scalar field, 4> m in, respectively as 

V eff {4>) = M 4 e ( ^ } + Ae LLK ^T-' p a- 3il+ ^ (2.17) 
,a+ 1 A 



mL„ = aA/ 4 +«(^ TT + — TT ) (2.18) 



and 



^ mm = I^(A[^(l + z)" 3 ( 1+ ^]^). 



(2.19) 



Note that 

m min is the inverse of the characteristic range of the chameleon force in a given medium. In the next 
section we best fit the model with the observational data. 

III. OBSERVATIONAL CONSTRAINTS 

We constrain the dimensionless coupling constant, /3, and model parameter a, in both exponential and power law 
potential, with the most recent observational data, SNe la, by employing the x 2 statistics. We assume that the 
energy density, p m , stands for the contribution from cold dark matter, 7 = 0. We also assume that M p i = 1, M = 1, 
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TABLE I: best fitted parameters 



Model 


a 


P 


x 2 


power law 


0.175 


4.51 


557.0747981 


exponential 


0.51 


2.31 


557.76250 




p p 

Fig. 1: The best-fitted confidence level for a and /3 for 
left)the power law potential and right) exponential potential 



A = 1. From Tabic U one finds the best fitted model parameters a and j3 in both cases of power law and exponential 
potentials. 

The confidence levels corresponding to the best fitted parameters are shown in FIG. 1). One sees that with 68.3%, 
95.4% and 99.7% confidence level the true values for both a ,(3 lie within the red, green and blue contours, respectively. 

Also, the best fitted distance modulus, p{z), with the experimental data in both exponential and power law potential 
cases illustrated in Fig. 2. 

With the best fitted model parameters, in the next section, we study in more details the chameleon mechanism by 
performing cosmological tests. 



IV. CHAMELEON MECHANISM AND COSMOLOGICAL TEST 

Complementary to the conventional chameleon mechanism given in pMj , here, we study the mechanism with respect 
to cosmological scales in different epoches. Just as an example, we take two values for the matter energy density, 
p m = I and p m = 3 where the corresponding size of the universe is a = ao = I (today) and a = 0.69 respectively. 
With the best fitted model parameters, in case of power law potential, the effective potential are shown in Fig. 3: 
top and bottom. From the graph, as matter energy density, p mi decreases or alternatively a — > ao, the minimum 
effective potential shifts to larger values of (j) m in- Moreover, the steepness of the effective potential, V e ff, near the 
minimum, also depends on p m ; a shallow minimum corresponds to a low chameleon mass. The mass of chameleon 
field, m, increases with p m . In other words, chameleon field has a greater contribution, as the universe approaches its 
earlier epoches. The same argument applies to the case of exponential potential case (Fig. 4). 

A cosmological quantity that usually considered as a significant parameter to measure the reliability and affectiveness 
of cosmological models is equation of state (EoS) parameter. In Fig. 5, a comparison between power law and 
exponential potential scenarios for the best fitted model parameters with observation is shown. The graph illustrates 
that while, up to about z ~ 1.5, the trajectories of effective EoS parameter in both scenarios overlap each other, 
they completely resolved afterward. Since there is no direct measurement of EoS parameter for the universe, one can 
not decide which form of potential is observationally privileged. The graph also shows that the current effective EoS 
parameter for the best fitted model parameters in both cases is u> e ff — —0.46. Next, wc examine our models with 
cosmological data for Hubble parameter. 

In Fig. 6 the best fitted Hubble parameter for both potential functions obtained from numerical calculation is 
compared with the observational data. As can be seen, for the redshift z < 0.5 the two curves overlap and relatively 
agree with the observational data. Within the distance 0.5 < z < 1.5, they again reasonably fit the data whereas have 
quite distinct shape. Nevertheless, similar to the case of EoS parameter, the two curves begin to resolve at z > 1.5. 
Again, the two models satisfy the observational data within the observationally fitted interval < z < 1.5. However, 
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at this stage, since the model is not tested against data for higher redshifts, z > 2, one still can not decide which of 
the two scenarios is observationally preferable. 

Next, we choose to compare the two models against observational data for velocity drift 
distributed over 1.5 < z < 5. With the velocity drift given by 



44[ where the data 



v = cH, 



o 



cH{z) 

TT7 : 



(4.1) 



a comparison of the two models with the observational data for high redshift are shown in Fig. 7. 

As can be seen from the graph one can choose power law potential function a more natural choice over exponential 
potential function. This is also in agreement with the discussion given in [45j . 



V. SUMMARY AND REMARKS 



This paper examines the chameleon mechanism in cosmological scales by taking into consideration power-law and 
exponential potential functions. The work is supported by best fitting the model parameters with the observational 
data. Departing from conventional approach to chameleon mechanism, we revisit the model of interacting chameleon 
field in the framework of cosmological scales. The result shows that in both cases of power-law and exponential 
potential the mass of interacting chameleon field reduces with the matter energy density. The two scenarios of power- 
law and exponential potential functions directly and indirectly are tested against observational data. We find that 
in the range of < z < 1.5, the dynamic of the effective EoS parameter is independent of the underlined models. 
However the circumstance is different in higher redshifts, i.e. z > 1.5. The EoS parameter also shows that the 
universe begins to accelerate at about z ~ 0.2 and currently is in quintessence era. We perform two observational 
tests comparing the dynamics of the model in two scenarios. First, the best fitted Hubble parameter derived from 
numerical computation in both cases are compared with the observational data. While both models are relativity in 
good agreement with the data for z > 1.5, no comparison can be made with the data for high redshifts. Second, we 
compare the best fitted velocity drift computed in both scenarios with the data for high redshift, 1.5 < z < 4.5. The 
test shows a better agreement between the data and best fitted velocity drift in power law potential case. 



[1] A.G. Reiss et al, Astron. J. 116, 1009 (1998) 



Fig. 4: The best fitted behavior of chameleon effective potential, V e ff((j>), in case of power law potential for (Top)p m = 1 and 

(Bottom)p m = 3. 



BcstfUofSNela 




Fig. 5: The best fitted effectivce EoS parameter for power law and exponential potential. 



[2] C. I. Bennet et al, Astrophys J. Suppl. 148:1, (2003) 

[3] A. G. Riess, et al, [Supernova Search TeamCollaboration] Astron J. 116 1009 (1998) 

[4] A. C. Pope, et. al, Astrophys J. 607 655, (2004) 

[5] D. N. Spergel, et. al., Astrophys J. Supp. 148 175, (2003) 

[6] V. Sahni, A. Starobinsky, Int. J. Mod. Phys. D 9 373-444, (2000) 

[7] S. Weinberg, Rev. Mod. Phys. 61 1(1989) 

[8] R. R. Caldwell, R. Dave and P. J. Steinhardt, Phys. Rev. Lett. 80 1582,(1998) 

[9] I. Zlatev, L. Wang and P. J. Steinhardt, Phys. Rev. Lett. 82 896,( 1999) 
[10] R. R. Caldwell, M. Kamionkowski, N. N. Weinberg, Phys. Rev. Lett. 91 071301, ( 2003) 
[11] Z. H. Zhu, M. K. Fujimoto and X. T. He, Astrophys J. 603 365-370, ( 2004) 
[12] Z. H. Zhu and J. S. Alcaniz, Astrophys J. 620 7-11 ( 2005 ) 

[13] Sadeghi, J., Setare, M. R., Banijamali, A., & Milani, F. 2008, Phys. Lett. B 662 92; Sadeghi, J., Setare, M. R., Banijamali, 
A., & Milani, F. 2009, Phy. Rev. D 79 123003; Guo, Z. K. et al. 2005, Phys. Lett. B 608, 177; Xia, J.-Q., Feng, B., & 
Zhang, X. 2005, Mod. Phys. Lett. A 20 2409; Setare, M. R. 2006, Phys. Lett. B 641 130; Zhao, W. & Zhang, Y. 2006, 
Phy. Rev. D 73 123509; Zhao, G.-B., Xia, J.-Q., Feng, B., & Zhang, X. 2007, Int. J. Mod. Phys. D 16, 1229; Setare, M. 
R., Sadeghi, J. & Amani, A.R. 2008, Phys. Lett. B 660 299; Setare, M. R. & Saridakis, E. N. 2008, Phys. Lett. B 668 177; 
Setare, M. R. & Saridakis, E. N. 2009, Int. J. Mod. Phys. D 18 549; Setare, M. R. & Saridakis, E. N. 2008, J. Cos. Astro. 
Phys. 09 026. 

[14] Cai, Y. F., Qiu, T., Piao, Y. S., Li, M. & Zhang, X. 2007, JHEP 0710 071. 

[15] Farajollahi, H., Salehi, A., Tayebi, F., Ravanpak, A. 2011, J. Cos. Astro. Phys. 05, 017. 

[16] Capozziello, S., Cardone, V. F., Carloni, S. & Troisi, A. 2006, Int.J.Mod.Phys. D15 69; 2003, Int.J.Mod.Phys. D12 1969. 
[17] M. R. Setare, Phys. Lett. B644:99-103,(2007) 

[18] Setare,M. R. & Jamil,M. 2010, Phys. Lett. B 690 1-4 ; Davis,A. C, Schelpe, C. A.O., Shaw, D. J., 2009, Phy. Rev. D 
80 064016 ; Ito, Y. & Nojiri, S. 2009, Phy. Rev. D 79:103008; Tamaki,T. & Tsujikawa,S. 2008, Phy. Rev. D 78 084028 ; 
Farajollahi, H. & Salehi, A. 2010b Int. J. Mod. Phys. D19:621-633; Mota,D.F. & Shaw, D.J. 2007, Phy. Rev. D 75, 063501; 
Dimopoulos, K. & Axenides, M. 2005, J. Cos. Astro. Phys. 0506:008. 

[19] D. H. Lyth and A. Riotto, Phys. Rept. 314, 1 (1999) 

[20] C. Wetterich, Nucl. Phys. B302, 668 (1988) 

[21] P. J. E. Peebles and B. Ratra, Ap. J. 325 ,L17 (1988) 

[22] R. R. Caldwell, R. Dave and P. J. Steinhardt, Phys. Rev. Lett. 80, 1582 (1998) 
[23] D. A. Easson, JCAP 070 2, 004 (2007) 

[24] T. Damour, G. W. Gibbons and C. Gundlach, Phys. Rev. Lett, 64, 123 (1990) 
[25] M. R. Setare, Elias C. Vagenas, Int. J. Mod. Phys. D18:147-157 (2009) 
[26] S. M. Carroll, Phys. Rev. Lett. 81 3067(1998) 

[27] S. M. Carroll, W. H. Press and E. L. Turner, Ann. Rev. Astron. Astrophys, 30, 499 (1992) 
[28] T. Biswas, R. Brandenberger, A. Mazumdar and T. Multamaki. Phys.Rev. D74 063501, (2006) 
[29] J. P. Uzan, Rev. Mod. Phys. 75, 403 (2003) 



7 



//' : 



Best fit of SNe la 




A 1 <W° 






v A* 




































1 J 1 4 1 



Fig. 6: The best fitted Hubble parameter for power law and exponential potential against observational data. 
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Fig. 7: The graph of best fitted velocity drift v for exponential and 
power law potentials in comparison with the data 



[30] B. Bertotti et al. Nature 425, 374 (2003) 

[31] G. F. Chew and S. C. Frautschi. Phys. Rev. Lett. 7, 394 (1961) 

[32] T. Damour, F. Piazza and G. Veneziano, Phys. Rev. D 66 , 046007 (2002) 

[33] S. Nojiri, S. D. Odintsov, Mod. Phys. Lett. A 19:1273-1280 (2004) 

[34] J. Khoury and A. Weltman: Phys. Rev. Lett. 93:171104,(2004) 

[35] D. F. Mota, J. D. Barrow, Phys. Lett. B581 141-146(2004); 

[36] J. Khoury and A. Weltman, Phys. Rev. Lett. 93,171104 (2004) 

[37] Ph. Brax, C. van de Bruck, A. C. Davis, J. Khoury and A. Weltman. Phys. Rev.D70, 123518 (2004) 
[38] C. Wetterich, Astron. Astrophys. 301, 321 (1995) 

[39] T. Damour and A.M. Polyakov, Nucl. Phys. B423, 532 (1994); Gen. Rel. Grav. 26, 1171 (1994) 
[40] G. Huey, P.J. Steinhardt, B. A. Ovrut and D. Waldram. Phys. Lett. B 476, 379 (2000) 
[41] C.T. Hill and G. C. Ross, Nucl. Phys. B311, 253 (1988) 

[42] J. Ellis, S. Kalara, K.A. Olive and C. Wetterich, Phys. Lett. B 228, 264 (1989) 
[43] D. F. Mota and C. van de Bruck, Astron. Astrophys. 421,71 (2004) 
[44] Liske J. et al, Mon. Not. Roy. Astron. Soc. 386, 1192 (2008) 

[45] J. Khoury, A. Weltman, Phys.Rev.D69:044026, (2004); H. Jason, et al., Phys.Rev.Lett. 105:261803,(2010); P.Brax, et al., 
Phys. Rev.D.82:083503, (2010); J. Khoury, A. Weltman, Phys. Rev. Lett. 93:171104,(2004) 



